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SECTION  I 


INTRODUCTION 

In  order  to  monitor  the  condition  of  each  aircraft  engine  in  its 
fleet,  the  Air  Force  periodically  samples  the  lubricating  oil  and  makes 
a  spectrometric  analysis  of  the  metal  content.  These  measurements  serve 
as  a  guide  to  identify  abnormal  wear  of  parts  and  to  schedule  inspections 
or  overhauls.  As  presently  operated,  the  Spectrometric  Oil  Analysis 
Program  (SOAP)  prescribes  alert  levels  for  the  concentration  of 
critical  metals  for  each  engine  type,  and  mandates  particular  actions  in 
the  field  when  such  levels  are  detected.  Pre-emptive  time-scheduled 
maintenance  is  performed  at  long  intervals  on  engines  not  flagged  in  the 
interim  by  SOAP. 

The  efficiency  of  the  SOAP  system  is  assessed  by  means  of  an  after- 
the-fact  classification  of  each  sample  analyzed.  Four  categories  are 
defined: 

a.  Failure.  An  oil-wetted  part  failed  although  the  immediately 
preceding  spectrometric  analysis  was  within  prescribed  limits. 

b.  Hit.  The  spectrometric  analysis  fell  outside  prescribed 
limits  and  the  immediately  subsequent  inspection  of  the 
engine  revealed  damage  to  a  part. 

c.  Routine.  The  spectrometric  analysis  fell  within  limits  and 

no  engine  malfunction  occurred  prior  to  the  next  spectrometric 
analysis . 

d.  Miss.  The  spectrometric  analysis  fell  outside  limits  but 
no  damage  was  found  upon  Inspection. 

Some  two  million  samples  are  analyzed  each  year  at  more  than  a  hundred 

laboratories  world-wide.  Of  these,  the  vast  bulk  are  routines;  all  other 

categories  are  comparatively  rare.  Failures  run  about  100-200  annually, 


hits  1200-1300,  and  misses  50-60.  Yet  the  SOAP  program  pays  for  its 
considerable  cost  (estimated  at  around  $10  per  analysis)  by  the  detection 
of  the  relatively  small  fraction  of  hits,  each  of  which  represents  the 
averting  of  a  possible  failure  that  could  result  in  the  loss  of  an  engine, 
an  airplane,  or  the  lives  of  aircrew.  Prevention  of  failure  is  the 
central  purpose  of  the  program,  and  the  goal  of  research  and  development 
in  SOAP  technology  is  to  decrease  the  ratio  of  failures  to  hits.  Even 
small  improvements  in  this  ratio  exert  a  large  favorable  leverage  on  the 
cost-effectiveness  of  SOAP. 

Various  supplementary  methods  of  analysis  are  under  investigation 
with  the  aim  of  characterizing  oil  samples  more  completely  to 
fine-tune  the  capability  of  SOAP  with  respect  to  identifying  incipient 
failure.  One  such  technique  is  particle  counting,  which  defines  the 
size  distribution  of  wear  metals  and  other  particulate  debris  if  present. 
This  report  develops  a  rationale  whereby  the  measured  total  volume  of 
particulates  in  a  sample  is  compared  with  an  expected  value  inferred  from 
spectrometric  analysis  of  the  same  sample.  Experiments  on  authentic  oil 
specimens,  using  this  method  of  screening,  indicate  that  the  effectiveness 
of  SOAP  in  preventing  engine  failure  can  be  increased  about  sixfold  by 
supplementing  each  spectrometric  analysis  with  a  particle  count. 


SECTION  II 


EXPERIMENTAL  PROCEDURE 


Equipment 

The  particle  counter  used  in  this  work  is  a  commercially  available 
light  scattering  laser-beam  instrument  (Spectrex  ILI-1000).  A  helium- 
neon  laser  beam  (wavelength  0.6328  microns)  scans  the  sample  in  a  circular 
path  for  a  measured  period  of  time,  and  a  photodetector  collects  pulses 
scattered  forward  by  suspended  particles.  The  intensity  of  each  pulse 
is  proportional  to  the  square  of  the  diameter  of  the  scatterer,  which  is  assumed 
spherical.  The  optics  are  designed  to  focus  only  pulses  which  originate 
within  a  zone  of  limited  length  (about  2  centimeters),  the  boundaries  of 
which  lie  within  the  sample.  This  length,  together  with  the  diameter  of 
the  beam  and  the  fixed  circular  rate  of  beam  travel,  determines  the 

O 

volumetric  scan  rate,  which  amounts  to  0.6165  cm  per  second  for  our 
particular  instrument.  The  counter  is  coupled  to  a  15-channel  pulse-height 
discriminator  which  sorts  and  stores  counts  according  to  their  intensity — 
thus,  according  to  the  size  of  the  scatterer  generating  the  pulse.  An 
additional  summation  channel  counts  all  pulses  received  by  the  15  sorting 
channels.  Particles  smaller  than  1-2  microns  do  not  scatter  light  at  the 
laser  wavelength  and  cannot  be  detected. 

In  order  to  avoid  multiple  scattering  of  pulses  and  also  keep  the 
rate  of  arrival  of  pulses  within  the  resolving  power  of  the  detector,  it  is 
required  that  the  total  concentration  of  particles  not  exceed  about  1500 

O 

per  cm  .  Further,  the  pulse  intensity  received  at  the  detector  from  a 
particle  of  a  given  size  is  dependent  upon  the  attenuation  (if  any)  of  the 
laser  beam  due  to  absorption  by  the  liquid  under  examination;  a  correction 
for  this  effect  can  be  made  if  the  opacity  of  the  specimen  is  independently 
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measured,  and  this  must  be  done  as  a  part  of  the  counting  experiment.  An 
opacity  of  35%  is  about  the  limit  for  convenient  operation;  above  this 
limit,  the  detector  rapidly  and  progressively  goes  blind,  except 
to  strong  pulses  (from  large  particles).  For  these  reasons  it  is  necessary 
to  dilute  oil  samples  before  counting,  since  both  the  particle  concentration 
and  the  opacity  of  most  used  engine  oils  are  well  beyond  the  bounds  mentioned. 

Our  particular  assembly  of  equipment  was  calibrated  before  use  against 
a  standard  polydisperse  suspension  of  polystyrene  spheres  in  water,  the 
total  concentration  and  size  distribution  of  which  were  traceable  to  NBS 
determinations.  Table  1  summarizes  the  results,  showing  the  boundaries 
(in  microns)  of  each  of  the  fifteen  channels  as  a  function  of  opacity. 

This  calibration  is  peculiar  to  the  particular  items  of  hardware  and  to 
the  specific  bias  applied  to  the  pulse-height  discriminator;  the  design 
of  the  manufacturer's  equipment  does  not  permit  any  parameter  of  the  system 
to  be  varied  independently  of  all  the  rest,  although  a  different  biassing 
would  result  in  a  different  calibration.  The  boundaries  shown  in  Table  1 
adequately  cover  the  range  of  particle  sizes  found  in  engine  oils. 

Technique 

A  procedure  was  developed  by  means  of  which  a  useful  particle  count 
could  be  accomplished  in  a  total  time  of  25  minutes  with  the  general-purpose 
counter  just  described.  Although  this  is  too  time-consuming  for  use  as  a 
field  method,  it  was  quite  satisfactory  for  the  intended  purpose  of  assessing 
the  usefulness  of  particle  counting.  Adaptation  to  field  use  appears  to  be 
a  straightforward  problem;  this  matter  will  be  discussed  later. 

In  the  present  work,  the  oil  sample  to  be  examined  was  diluted 
to  reduce  the  opacity  and  particle  concentration  to  tolerable  values.  As  a 
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TABLE  1.  CALIBRATION  OF  PARTICLE  COUNTER 


Spectrex  I L I - 1 OOO  Serial  No.  131,  Coupled  to  Particle 
Profile  Attachment  Serial  No.  0038 

Bias  Settings:  Gain  =  500,  Threshold  =  260 


Entries  show  particle  diameters,  in  microns,  which  bound 
each  channel  as  a  function  of  opacity 


[ 

1 

Channel 

Opaci ty 

7.5 

10 

15 

20 

25 

30 

35 

2.07 

2.16 

2.34 

2.53 

2.74 

2.96 

3.21 

01 

5.22 

5.42 

5.86 

6.32 

6.83 

7.37 

7.96 

02 

10.67 

11.07 

11.94 

12.86 

13.87 

14.94 

16.11  : 

03 

16.17 

16.77 

18.06 

19.44 

20.93 

22.53 

24.26 

04 

24.36 

25.26 

27.16 

29.20 

31 .41 

33.78 

36.32 

,  05 

1 

32.07 

33.24 

35.71 

38.37 

41  .25 

44.32 

47.62 

!  06 

36.26 

37.57 

40.36 

43.35 

46.58 

50.04 

53.75 

i  07 

i 

39.82 

41  .26 

44.31 

47.58 

51  .11 

54.89 

58.94 

!  08 

43.01 

44.56 

47.84 

51  .36 

55.17 

59.23 

63.59 

09 

■ 

45.88 

47.52 

51  .02 

54.76 

58.81 

63.14 

67.77 

10 

48.72 

50.46 

54.16 

58.13 

62.42 

66.99 

71.90 

11 

51  .37 

53.20 

57.10 

61.27 

65.78 

70.59 

75.75 

i  12 

i 

54.15 

56.08 

60.17 

64.56 

69.30 

74.37 

79.79 

13 

56.68 

58.69 

62.97 

67.56 

72.51 

77.80 

83.46 

14 

59.45 

61  .56 

66.04 

70.84 

7C.02 

81  .55 

87.48 

15 

62.27 

64.47 

69.16 

74.17 

79.59 

85.37 

91  .57 

o 

diluent,  we  used  a  fixed  volume — usually  12.5  cm — of  heavy  white  paraffin 
oil  (opacity  7.5%,  due  entirely  to  reflection  at  the  cell  windows)  contained 
in  a  conventional  spectrophotometer  cell  of  5  cm  path  length.  The  engine 
oil  was  pipetted  in  to  produce  a  suspension  of  known  dilution,  usually 
about  70  fold;  in  exceptional  cases  a  lower  dilution  may  be  used  to  accomodate 
an  oil  having  very  low  particle  concentration.  The  mixture  was  well  shaken 
and  then  degassed  for  a  minute  or  two  at  low  (house)  vacuum  in  a  desiccator. 
The  opacity  was  then  measured  with  an  accessory  photometer.  A  single  run 
consisted  of  a  preliminary  blank  count  on  the  paraffin,  followed  by  a  second 
count  after  addition  of  the  engine  oil.  The  counting  time  was  450  seconds 
in  each  case — found  by  experience  to  be  about  the  minimum  required  to  give 
reproducible  counting  statistics  in  the  highest  populated  channels.  The 
viscosity  of  the  paraffin  (about  80  centipoise)  was  sufficient  to  obviate 
any  appreciable  sedimentation  during  the  counting  period.  Counts  stored 
in  the  pulse-height  discriminator  for  both  blank  and  sample  were  printed 
out  in  hard  copy  for  each  of  the  channels;  these  data,  along  with  the 
measured  dilution  and  opacity,  were  fed  off-line  to  a  computer  and  worked 
up  by  a  program  called  NEWPART. 

The  NEMPART  program 

This  computer  program  has  been  developed  for  the  purpose  of  transforming 
the  raw  particle-count  data  into  a  table  showing,  for  round  values  of  the 
particle  diameter  at  5-micron  intervals,  two  quantities  N  and  P,  defined 
as  follows; 

3 

N  =  number  of  particles,  per  cm  of  oil  as  received, 
which  are  larger  than  the  listed  diameter 

P  =  volume  of  particles,  microliters  per  liter  of  oil 

as  received,  which  are  larger  than  the  listed  diameter 
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Figure  2.  Newpart 


Certain  mathematical  constants  generated  during  the  data  reduction  are  also 
recorded.  Figures  1  through  3  are  typical  NEWPART  printouts  for  authentic  oil 

Table  2  gives  a  quick  overview  of  the  functions  of  this  program,  and 
Appendix  A  Is  a  complete  listing  In  FORTRAN  language.  A  few  comments  may 
be  helpful  in  illuminating  essential  details. 

The  Input  is  supplied  on  punched  cards.  A  parity  check  assures  that 
the  sum  of  all  channel  counts  equals  the  total  number  of  counts;  this 
provides  a  safeguard  against  errors  in  keypunching. 

The  data  for  blank  and  sample  are  separately  smoothed  by  a  logarithmic 
fit  of  the  function. 

N  »  Noexp(-aDk)  +  N0exp(-aDi5k)  (1) 

which  contains  three  adjustable  constants  (NQ,  a,  and  k)  and  is  formally 
similar  to  a  Wei  bull  distribution.  However,  we  must  regard  it  only  as  an 
empirical  equation  which  fits  available  data  well ---a  graphic  example  is 
shown  in  Figure  6  later  on — and  marshals  them  in  a  convenient  form  for 
analysis.  The  second  term  on  the  right  arises  from  the  fact  that  the  top 
of  channel  15  is  closed,  rather  than  open  ended,  so  that  particles  larger 
than  this  top  limit  will  not  be  counted  even  if  present;  but  they  will 
presumably  lie  on  the  curve. 

With  respect  to  the  low  cutoff  diameter  Dj_q ,  we  take  this  either  as  one 
micron  (the  smallest  dimension  visible  to  the  particle  counter)  or  as  the 
smallest  integer  diameter  for  which  the  difference  between  the  sample  and 
blank  curves  is  positive,  whichever  is  greater.  Extrapolation  of  the  fitted 
curves  to  regions  of  smaller  D  generally  yields  specious  values  (predicts  too 
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high  a  concentration);  in  particular,  N0  is  not  a  physically  meaningful 
estimate  of  the  total  number  of  particles. 

The  high  diameter  cutoff  D^j  expresses  the  fact  that  no  real  sample 
contains  indefinitely  large  particles.  Counts  in  the  highest  channels  are 
rare  events,  governed  by  Poisson  statistics  in  which  the  standard  errors 
of  Ng  and  N$  approach  the  square  roots  of  the  respective  quantities.  We 
may  then  apply  the  conventional  criterion 

z  »  (ns-nb)/(ns+nb)1/2  (2) 

to  determine  whether  or  not  a  significant  difference  exists  between  sample 
and  blank.  We  choose  (arbi trari ly )  to  declare  the  difference  in  concentration 
as  insignificant  at  the  50%  level,  for  which  Z  =  0.675.  A  different  criterion 
would  raise  or  lower  DHj  slightly,  but  the  choice  is  not  critical. 

Once  Dlo  and  DHI  have  been  determined,  the  quantity  P(D)  is  computed 
for  round  values  of  D  at  5-  micron  intervals  in  the  range  Dlo<d<dHI  as 
indicated  in  the  synopsis.  The  integration  may  be  made  by  straightforward 
numerical  quadrature  or  by  noting  that  each  of  the  two  integrals  (blank  and 
sample)  may  be  expressed  in  terms  of  incomplete  gamma  functions: 

D3(dN/dD)dD  =  -  N„  a'3/k  I"  r  (1+3/k)  -  T  (l+3/k)l  (3) 

LaDk  aDf0  J 


Efficient  algorithms  are  available  to  evaluate  the  gamma  functions  and  these 
are  used  by  NEWPART  rather  than  the  more  cumbersome  numerical  integration. 
The  additional  term  in  N^j  Djj  is  simply  a  rough  estimate  of  the  volume 
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(usually  small)  of  those  particles  which  may  be  distributed  in  some  undefined 
way  above  the  high  cutoff  diameter;  there  is  a  50%  probability  that  this 
increment  is  genuine  for  any  given  sample. 


It  is  well  to  emphasize  that  NEWPART  as  listed  contains  the  instrument 
calibration  data  shown  in  Table  1.  If  a  different  instrument  were  used,  or 
if  the  bias  were  changed  on  the  same  instrument,  the  appropriate  data  would 
have  to  be  introduced  into  NEWPART  to  replace  the  present  tabulation. 

Finally,  note  Figure  3,  which  states  NO  SIGNIFICANT  PARTICLE  COUNT  IN 
THIS  SAMPLE.  The  computer  program  returns  this  message  in  any  case  for 
which  the  fitted  blank  curve  (N  vs  D)  lies  everywhere  above  the  sample 
curve,  or  for  which  fy|i  <  5  microns  with  a  correspondingly  small  and 
uncertain  maximum  value  of  P  corresponding  to  —  usually  less  than  .05. 
This  nil  finding  occurs  rather  frequently  in  the  examination  of  authentic 
engine  oils. 

Synthetic  suspensions 

In  order  to  test  the  validity  of  both  the  experimental  techniques  and 
the  data  reduction  program,  measurements  were  made  on  synthetic  suspensions 
having  known  properties. 

One  such  experiment  was  to  make  a  normal  450-second  count  on  the 
standard  suspension  of  polystyrene  spheres  against  which  the  instrument  was 
originally  calibrated.  This  was  fundamentally  a  test  of  the  computer  program, 
although  it  also  would  reveal  any  outright  blunders  in  performing  the 
calibration.  Figure  4  Is  a  NEWPART  printout,  using  an  artificial  blank  <nput 
of  negligible  particle  content.  Superimposed  on  the  printout  is  a  table  of 
the  known  values  of  N  vs  0  for  the  standard,  together  with  values  of  P 
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derived  by  numerical  Integration.  These  data  are  seen  to  agree  well  with 
the  corresponding  values  generated  by  NEWPART.  The  concordance  could 
presumably  be  improved  further  by  increasing  the  counting  time,  but  this 
test  Indicates  that  the  normal  450-second  run  is  adequate  for  practical 
purposes . 

A  second  kind  of  test  consists  of  counting  a  suspension  for  which  the 
total  volume  of  particles  is  directly  known — prepared  by  weighing  out  a 
quantity  of  pure  metal  powder  and  dispersing  it  into  a  known  mass  (or 
volume)  of  clean  paraffin  oil.  A  specimen  of  this  suspension,  diluted 
if  necessary,  is  counted;  the  value  of  P  (for  DLq)  derived  from  the  particle 
count  is  then  compared  with  the  value  which  is  obtained  by  a  trivial 
calculation  involving  the  known  densities  of  the  paraffin  and  the  metal  in 
question.  This  experiment  has  been  performed  many  times  with  uniformly 
satisfactory  results;  the  NEWPART  estimate  consistently  lies  within  +3% 
of  the  true  value.  Table  3  shows  the  results  of  one  such  experiment  which 
is  particularly  interesting  because  it  was  a  blind  test;  the  suspension 
was  prepared  by  a  colleague,  and  neither  the  mass  concentration  nor  the 
identity  of  the  metal  was  revealed  until  after  the  particle  count  (Figures  ) 
had  been  made.  Figure  6  is  a  plot  of  the  raw  values  of  N  for  sample  and 
blank  showing  the  best- fitted  particle  size  distribution  curves  (Equation  1) 
and  the  cutoff  diameters  generated  in  the  NEWPART  analysis. 
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TEST  ON  SYNTHETIC  SAMPLE 


synthetic  sample**** 
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Figure  5.  Analysis  o.f  Synthetic  Suspension 
of  magnesium  powder  in  paraffin  oil 
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Figure  6.  Raw  Data  for  Synthetic  Suspension 
of  magnesium  powder  In  paraffin  oil,  showing 
best- fitted  curves  generated  by  NEWPART  program 
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SECTION  III 


EXAMINATION  OF  AUTHENTIC  OILS 


Data  base 

Over  the  past  18  months,  AFML  has  received  periodic  inputs  of  authentic 
oil  samples  collected  in  the  field,  together  with  the  spectrometric  metal 
analysis  made  by  a  SOAP  laboratory.  These  samples  include  all  SOAP  failures 
and  hits  experienced  at  the  several  participating  Air  Force  stations,  plus 
a  random  selection  of  routines,  comprising  in  all  a  stock  of  about  1000 
specimens  of  oil.  Out  of  this  grand  lot,  a  sub-lot  of  200  samples  was 
subjected  to  particle  counting  by  the  procedure  already  described. 

The  essential  data  from  these  experiments  are  collected  for  reference 
in  Appendix  B.  All  available  failures  (31)  were  examined,  the  remainder 
being  hits  (24)  and  routines  (145)  randomly  chosen;  no  misses  were  received 
from  the  field.  The  samples  analyzed  represent  a  variety  of  aircraft  and 
engine  types;  the  service  age  of  the  oils  at  the  time  of  sampling  covers 
a  wide  range. 

For  each  sample,  we  list  two  parameters,  X  and  Y,  which  form  the  basis 
of  the  screening  system  to  be  described  presently.  X  is  the  total  value 
of  P,  corresponding  to  Dj_q ,  as  reported  by  the  NEWPART  program.  Y  is  an 
independent  estimate  of  P  inferred  from  the  spectrometric  analysis  of  the 
same  sample,  given  by 

Y  =  p^miVi  (4) 

where  m.  and  V.  are  respectively  the  spectrometric  mass  concentration  (ppm) 
and  the  specific  volume  of  the  metal  i,  and  p  is  the  density  of  the  lubricating 
oil  (0.85  gm/cm^  to  a  sufficiently  good  approximation).  The  following  tabulation 
shows  the  specific  volumes  (cm^/gm)  of  the  twelve  metals  customarily  determined 
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in  SOAP  analyses. 


A1 

0, 

,37 

Pb 

0.09 

Si 

0.43 

Cr 

0, 

,14 

Mg 

0.57 

Ag 

0.10 

Cu 

0. 

,11 

Mo 

0.10 

Sn 

0.14 

Fe 

0. 

,13 

Ni 

0.11 

Ti 

0.22 

Metal  particles  in  oil  are  of  course  usually  present  as  alloys  rather  than 
in  pure  form.  Equation  4  contains  the  implicit  assumption  that  the  specific 
volume  of  an  alloy  is  a  linear  function  of  this  property  for  its  constituent 
metals . 

To  make  this  matter  clear  by  a  specific  example,  refer  to  the  entry  for 
sample  F-035  in  Appendix  B  and  to  Figure  1,  which  is  the  particle  count  for 
this  sample.  The  entry  for  X  in  Appendix  B  is  6.83,  which  is  simply  the 
value  of  P  read  from  Figure  1  at  D^q.  To  compute  Y,  we  must  know  the 
spectrometric  analysis  for  this  sample.  The  result  reported  by  the  SOAP 
laboratory  (ppm  by  mass)  was  A1  =  5 ,  Cr  =  1 ,  Mg  =  1 ,  and  Si  =  10,  no  other 
metals  being  detected.  Applying  Equation  4  gives  Y  =  5.83  as  shown  in 
Appendix  B.  Entries  for  all  other  samples  are  the  result  of  the  same  process 
requiring  both  a  particle  count  and  a  spectrometric  analysis. 

Interpretation  of  Data 

We  now  construct  a  rationale  to  show  how  particle-count  data  can  be 
applied  to  supplement  spectrometric  analysis  and  improve  SOAP  decision-making 
The  scheme  rests  on  the  comparison  of  X  and  Y  values  for  the  samples  in 
question.  If  the  ratio  X/Y  exceeds  unity,  the  sample  is  considered  defective 
whether  or  not  established  SOAP  protocols  mandate  an  alarm.  The  argument 


A  value  X/Y>1  constitutes  evidence  for  the  existence  of  one  or  both  of 


two  conditions: 


a.  Large  particles  of  metal  are  being  generated,  above 
10-15  microns  in  diameter,  which  is  the  upper  limit 
of  effective  spectrometric  analyzability  for  most 
metals.  Appearance  of  such  large  metal  particles, 
undetected  by  spectrometry ,  is  generally  indicative 
of  the  onset  of  shaving,  milling,  or  fatigue 
phenomena  which  are  not  normal  wear  modes. 

b.  The  oil  contains  substantial  quantities  of  non-metallic 
debris  which  in  no  event  could  be  detected  by  spectro¬ 
metric  analyses.  Such  particles  may  originate  from 
degradation  of  the  oil  itself  through  polymerization 

or  sludge  formation,  or  from  external  contamination  by 
dirt,  water,  or  damage  to  non-metallic  seals  and  gaskets. 
Presence  of  such  debris  is  an  early  warning  of  sub¬ 
standard  oil  condition  presaging  eventual  increases  in 
wear-metal  production. 


On  the  other  hand,  if  X/Yd  and  if  the  spectrometric  analysis  lies  within 
established  limits,  normal  wear  is  indicated.  The  particle  size  distribution 
of  wear  metal  (and  of  non-metallic  debris,  if  present)  is  such  that  most 
particles  do  not  exceed  1-2  microns  in  linear  dimension — the  smallest  size 
of  object  visible  to  the  particle  counter.  Therefore,  oil  and  engine  are 
operating  in  a  healthy  condition. 

In  brief,  it  is  suggested  that  the  existing  single  spectrometric  test 
should  be  augmented  by  the  introduction  of  a  second  screen  (particle  counting). 
The  oil  sample  is  then  required  to  pass  both  tests.  Figure  7  shows  the  scheme 
graphical ly. 

Results  and  discussion 

The  feasibility  of  the  dual -screening  principle  is  readily  demonstrated. 
Figure  8  is  an  X,Y  plot  of  the  extant  data,  showing  that  spectrometric  routines 
(open  circles)  fall  predominantly  below  the  curve  X  =  Y,  whereas  most  failures 
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SPECTRO 


Figure  7,  Huai  Screening  Scheme 


2  •  )  LIE  IN  THIS 
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Y,  VALUE  OF  P 
FROM  SPECTROSCOPY 


Figure  8.  Failure  Prediction 
based  on  comparison  of  two  estimates  for  P, 
microliters  of  particulates  per  liter  of  oil 
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(filled  circles)  lie  above  it;  hits  (crossed  circles)  are  almost  evenly 
distributed.  Table  4  summarizes  the  results  in  a  compact  form. 

According  to  this  evidence,  supplementary  particle  counting,  had  it 
been  a  part  of  the  SOAP  screening  process,  would  have  identified  84%  of 
the  failure  samples  as  defective  and  would  have  caused  an  inspection  to  be 
made  before  the  failure  occurred.  In  the  case  of  hits,  almost  half  of  the 
samples  had  already  surpassed  the  X/Y  limit  by  the  time  the  spectrometri c 
analysis  gave  an  alarm;  the  source  engines  would  have  been  flagged  for 
inspection  at  some  previous  juncture.  And  if  the  spectrometric  routine  samples 
are  indeed  representative  of  their  class,  then  11%  of  all  engines  now  in 
service  contain  oil  of  a  suspicious  character;  these  are  the  engines  most 
likely  to  yield  future  spectroscopic  hits  or  failures.  Again,  supplementary 
particle  counting  would  have  precipitated  an  inspection  of  each  of  these 
engines  at  some  earlier  time,  resulting  in  either  a  hit  or  miss  finding.  A 
point  to  be  emphasized  is  that  this  would  have  occurred  the  first  time  the 
X/Y  ratio  exceeded  unity.  There  could  be  no  accumulation  of  dirty  oil  in  the 
fleet;  it  would  be  purged  upon  discovery. 

Ultimately,  i  f  we  are  to  make  reliable  cost  effectiveness  studies,  we 
need  to  know  how  supplementary  particle  counting  would  affect  both  the  annual 
number  of  failures  and  the  annual  number  of  inspections  (hits  plus  misses).  The 
actual  magnitude  of  these  parameters  can  be  determined  only  by  a  program  of 
real-time  dual  analysis.  Our  present  experiments  allow  us  to  judge  with  fair 
confidence  that  the  number  of  failures  will  be  reduced  to  about  16%  of  the 
present  level;  on-line  experiments  will  serve  principally  to  enlarge  the 
statistical  data  base  and  refine  the  estimate  a  few  percent  either  way.  The 
situation  is  quite  different  with  respect  to  forecasting  the  number  of  inspections 


25 


and  their  hit/miss  ratio.  Our  present  experiments  on  one-time  grab  samples 
can  afford  no  quantitative  answers  to  this  question,  which  is  of  critical 
importance . 

One  might  be  tempted  to  suppose,  from  the  results  shown  in  Table  4 
that  supplemental  particle  counting  would  single  out  about  11%  of  all 
spectroscopic  routine  samples  for  inspection — an  intolerable  burden  of 
cost  and  labor  even  if  failures  were  entirely  eliminated.  However,  this 
is  an  erroneous  conclusion.  All  we  can  say  is  that  11%  of  all  oil  now  in 
service  is  X/Y  defective;  we  cannot  connect  this  figure  with  an  anticipated 
on-line  inspection  rate  except  to  say  that  it  would  be  an  extravagant 
estimate.  Let  us  examine  the  matter  a  little  further  to  see  why  this  is  so. 

Referring  again  to  Appendix  B,  consider  the  data  showing  T,  the 
operating  time  since  oil  change,  for  spectroscopic  routine  samples.  Each 
of  these  samples  was  found,  by  particle  counting,  to  be  either  defective 
( X/ Y >1 )  or  non-defective  (X/Y<1).  If  we  compute  the  mean  and  the  variance 
of  T  for  each  of  these  subsets,  we  find  by  application  of  the  standard  t- 
and  F-  tests  that  there  is  no  reasonable  probability  that  they  form  two 
different  populations.  This  result  simply  confirms  that  all  the  routine 
samples  were  indeed  randomly  selected — i.e.,  the  particle  counts  were  made 
at  random  times  since  the  last  oil  change;  on  the  average,  the  defectives 
had  not  been  in  service  longer  than  the  non-defectives,  or  vice  versa.  Note 
now  that  every  engine  that  is  yielding  routine  samples  must  eventually, 
upon  continued  sampling,  issue  a  hit  or-  miss  sample,  unless  it  first  fails 
or  is  removed  from  service  for  time-scheduled  maintenance.  This  is  true 
whether  or  not  the  spectroscopic  analysis  is  supplemented  by  particle  counting. 
But  in  the  experiments  reported  here,  each  defective  sample — although  the 


source  engine  had  not  yet  failed  nor  generated  a  SOAP  alarm — came  from  oil 
that  had  been  defective  for  some  period  of  time.  How  long?  Did  damage 
already  exist?  What  happened  later  on?  We  cannot  say.  On  the  other  hand, 
each  non-defective  sample  originated  from  an  oil  stock  in  an  engine  which 
would  continue  to  pass  subsequent  X/Y  screening,  were  it  to  be  performed, 
for  some  additional  period  of  sampling.  How  long?  Again,  one-time  analysis 
affords  no  answer. 

Our  basic  premise  is  that  an  X/Y-defecti ve  sample  is  in  fact  a  genuine 
indicator  of  abnormal  engine  wear.  To  the  extent  that  this  concept  is  valid 
and  uncomplicated  by  extraneous  factors,  real-time  supplemental  particle 
analysis  would  provide  only  hit  findings,  resulting  in  early  detection  and 
correction  of  the  severe  wear;  the  engine  would  then  revert  to  normal  slow- 
wear  status.  Thereafter  it  might  pass  dozens  or  hundreds  of  dual  analyses. 
From  this  standpoint,  it  is  clear  that  the  size  of  the  X/Y-defecti ve  oil 
pool  now  in  the  fleet  is  no  index  for  estimating  a  new  on-line  inspection 
rate.  That  is  quite  another  thing:  namely,  the  number  of  inspections 
per  analysis.  If  all  X/Y-tri ggered  inspections  were  indeed  hits,  then  the 
increases  in  the  annual  number  of  inspections  would  approximate  closely  to 
the  decrease  in  the  number  of  failures  {roughly  a  hundred,  calculated  as 
84%  of  the  100-200  that  now  occur  in  a  year);  and  every  one  of  these 
additional  inspections  would  represent  pure  benefit.  We  recognize,  of  course, 
that  this  ideal  will  not  be  realized  in  practice.  Some  fraction  of  X/Y- 
mandated  inspections  on  spectroscopic  routines  will  be  misses;  these  false 
alarms  will  result  at  best  in  a  prophylactic  oil  change,  and  we  need  to 
know  the  magnitude  of  this  burden.  Only  continuous  on-line  supplemental 
analysis,  with  inspection  feedback,  can  provide  the  required  information;  it 
cannot  be  extracted  from  our  present  data. 
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Finally,  note  that  the  X/Y  criterion  Is  potentially  a  flexible  one 
insofar  as  the  critical  screening  value  is  concerned.  The  value  of  unity 
used  in  the  present  exploratory  experiments  is  not  necessarily  a  universal 
constant,  or  even  an  optimum  value.  Although  It  works  reasonably  well  when 
applied  to  the  limited  data  available,  we  must  remember  that  these  data — in 
addition  to  their  shortcomings  as  single  grab-samples — represent  an 
indiscriminately  lumped  pool  of  samples  from  different  types  of  aircraft 
and  engines;  these  in  turn  differ  greatly  among  themselves  with  respect  to 
severity  of  usage  and  rate  of  wear.  For  example,  fighter  engines  frequently 
operate  at  maximum  power,  while  transport  engines  are  seldom  run  at  maximum 
output.  Further,  depending  upon  the  location  of  the  engine  airscoop  relative 
to  ground  level,  some  aircraft  are  more  prone  than  others  to  ingest  dirt 
and  debris  which  may  find  its  way  into  the  oil  system.  For  these  reasons, 
it  is  reasonable  to  expect  that  in  future,  given  sufficient  oil  test  data 
identified  by  aircraft  type,  specific  X/Y  values  may  be  established  which 
optimize  screening  for  each  particular  case — less  than  unity  for  some  types, 
greater  for  others. 


Practical  counting  technique:  a  prospectus 

We  now  consider  briefly  the  problem  of  streamlining  the  particle  counting 
procedure.  As  a  prelude  to  any  immediate  program  of  real-time  testing,  it  is 
necessary  to  reduce  the  time  required  for  a  particle  count  to  the  same  order 
of  magnitude  (5  minutes)  as  that  for  a  spectroscopic  analysis.  This 
objective  can  be  achieved  through  straightforward  adaptions 
of  existing  equipment.  First,  the  particle  counter  must  be  modified  to 
accept  an  undiluted  thin  film  of  sample,  introduced  into  throwaway  cuvettes 
or  into  a  permanently  installed  flowthrough  cell  which  can  be  flushed  clean 
with  a  volatile  solvent  between  runs.  This  obviates  the  necessity  for  a 


separate  blank  count  and  eliminates  the  time  and  labor  associated  with  this 
extra  procedure.  The  film  must  be  sufficiently  thin  to  preclude  multiple 
scattering  and  to  produce  a  useable  opacity  of  35%  or  less.  Simple 
calculations  indicate  that  a  film  thickness  of  0.75  mm  will  guarantee  both 
these  conditions;  further,  the  total  counting  time  is  reduced  to  about  3 
minutes.  Refer  to  Appendix  C  for  details.  Preliminary  tests  on  neat  oils 
confirm  these  predictions  but  also  show  that  the  boundary  diameters  of  each 
channel  in  the  pulse-height  discriminator  are  quite  sensitive  to  the  exact 
location  of  the  thin  sample  within  the  2-cm  sensitive  zone  of  the  counter. 

The  counter  must  therefore  be  calibrated  for  a  particular  position  of  the 
sample  cell,  which  must  in  turn  be  rigidly  controlled  by  a  suitable  mechanical 
mounting . 

The  second  adaptation  is  to  provide  on-line  real-time  data  readout  by 
a  dedicated  microcomputer.  The  computer  program  would  be  a  simplified 
version  of  NEWPART,  purged  of  the  unnecessary  elements  connected  with  the 
manipulation  of  data  originating  from  the  blank  count;  additionally,  it 
might  be  augmented  to  accept  the  spectroscopic  analysis  as  an  auxiliary 
input  and  compute  X/Y  directly  as  a  sole  output,  with  the  full  NEWPART 
output  also  available  at  the  operator's  demand.  One  of  the  more  elaborate 
hand-held  programmable  calculators  with  thermal  printout,  interfaced  to 
the  pulse-height  discriminator ,  would  be  a  satisfactory  expedient. 

For  the  purpose  of  immediate  future  tests,  these  fundamental 
modifications  can  be  patched  onto  the  available  instruments  in  the 
laboratory.  In  the  longer  run,  such  features  would  best  be  integrated  into 
the  equipment  by  the  manufacturer,  along  with  additional  automation  to 
achieve  maximum  convenience  and  throughput.  A  few  examples  include 
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o  Built-in  measurement  of  opacity,  and  adjustment 
to  a  standard  value  (e.g.,  35X)  by  a  servo- 
controlled  variable  density  filter  rotated  into 
the  laser  beam. 

o  Open-ended  channels  in  the  pulse-height  discriminator, 
sorting  pulse  heights  directly  on  an  oversize  basis; 
independent  control  of  the  threshold  intensity 
(particle  diameter)  for  each  channel. 

o  Automatic  flush  and  feed,  to  allow  a  batch  of  samples 
to  be  run  serially  without  intervention  of  the  operator. 


Improvements  of  this  kind  represent  state-of-the-art  technology.  While 
custom  design  would  be  premature  at  the  present  stage  of  development,  it 
offers  capabilities  susceptible  to  quick  utilization  when  and  if  particle 
counting  is  adopted  as  a  supplement  to  SOAP. 
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SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  exploratory  phase  of  this  work  has  been  completed.  Particle 
counting  appears  to  be  a  promising  candidate  method  of  supplemental 
analysis  from  the  standpoint  of  failure  prevention.  The  question  of  cost 
effectiveness  remains  to  be  investigated,  and  this  calls  for  a  regime  of 
testing  different  from  what  has  been  conducted  thus  far. 

It  is  recommended  that  further  work  proceed  along  the  following  lines. 

a.  Modify  equipment  and  procedure  to  accomplish  rapid  particle 
counting  in  the  thin-film  mode,  incorporating  data  readout  from  a  dedicated 
microcomputer. 

b.  Using  this  nimble  capability,  established  in  a  central  R&D 
laboratory,  conduct  a  program  of  real  time  particle  analysis  on  every 
spectrometri c  sample  from  a  number  of  operational  aircraft  (perhaps  20). 

These  aircraft  should  be  selected  to  provide  a  cross  section  of  engine 
types  and  severity  of  service;  the  test  period  should  be  of  sufficient 
duration  (at  least  a  year)  to  guarantee  a  variety  of  maintenance  experience. 

c.  During  the  test  period,  issue  inspection  advice  and  receive 
inspection  feedback  in  order  to  assess  the  reliability  and  cost  effectiveness 
of  supplemental  particle  counting.  Continuously  refine  decision  criteria 
insofar  as  possible  to  improve  performance. 

d.  Assuming  favorable  findings,  move  wel 1 -automated  equipment 


and  necessary  instructions  out  to  one  or  more  SOAP  laboratories  to 
demonstrate  the  system  as  a  field  method. 


APPENDIX  A 


NEWPART  COMPUTER  PROGRAM 

The  following  pages  comprise  a  FORTRAN  listing  of  the  program.  It  is 
reproduced  here  primarily  for  the  reader  who  may  wish  more  detailed  information 
than  that  provided  in  the  synopsis  (Table  2  ).  For  the  reader  who  may  wish 
to  adapt  this  program,  or  parts  of  it,  for  use  in  connection  with  his  own 
particle-counting  equipment,  some  cautions  are  in  order.  First,  his  computer 
must  have  access  to  the  I SML  (International  Statistical  and  Mathematical 
Library)  in  order  to  call  out  algorithms  for  the  incomplete  gamma  function; 
lacking  this  capability,  he  must  modify  the  program  so  as  to  perform  the 
relevant  integrations  by  numerical  quadrature.  Second,  in  every  case  the 
BLOCK  DATA  section  must  be  replaced  with  information  pertinent  to  his  own 
equipment;  and  the  same  remark  applies  at  other  places  where  specific 
instrumental  constants  appear---for  example,  in  SUBROUTINE  OILSUB,  line  48, 
the  constant  1.622  is  the  time  required  for  our  instrument  to  scan  a  volume 
of  one  cubic  centimeter. 


PROGRAM  OIL 


1  PROGRAM  OILIIN PUT, OUTPUT, T  A*c  5  =  INPUT , T  APE6*0 DTP JT I 

DIMENSION  TITLE<12)  ,IV(1») 

COMMON/1  NPUT/UIi.,N,T,  OMEGA,  2(10,16)  ,I8LANK,PREC 
CO MM ON/OUTPUT/ D£lM(21> , DE L°(21),A2,FK2,FN[j, IT ER.IM, 10(211, MI 
S  COHMON/SOlVB/A8,F<3,FNOO,Z,  IDLOW, IOHIGH,KRIT 

CQHMON/GRI08/IT£RA,NO£N3  """"  '  ~  " — - 

NAMELIST/DATA/DIL,N,T,OMSGA,IBLANK,PREC,IZTAI 

1  FORMAT (12A6) 

2  FORMAT (16F5. 0)  '  " 

10  3  FORMAT </3X,»CHEC<  THIS  CART* FI  OK, 161b) 

4  FORMAT  (//12X,*SJMHARV  OF  RES  JL  TS»  F /9X  ,  *0*  ,  7X ,  »N* ,  1 1  < ,  »  P»/) 

5  ~ FORMAT (lib, 1P2E12. 2) 

6  FORMAT  (//4X,»N=NUMB£R  OF  PARTICLES, PER  CM3  0-  SAM»LE, WHICH  HAVE  A 

_ ZDIAMtTER  LARGER  THAN  0  */  LX, * P=20LUME  OF  PARTICLES, IN  PPM  OF  SAMP^ 

15  2  L  E  ,  WHICH  HAVE  A  DIAMETER  -ARGER  THAN  D* /4X ,  *0  =  3 t A 1£ T £R  0-  PARTIC 

2LES  IN  MICRONS*) 

7  FORMAT  </3x,*A  OR  <  IS  EDUAl  OR  LESS  THEN  ZERO.*) 

9  FORMAT (/3X,*LESS  THEN  THREE  DATA  POINTS,*) 

10  F0RMAT(1X,*A(8LA2<  ) =*,1»E9. 3, 3X,*K(BLAN<  ) =  * , E9.  3 , 5 < , » ND * » , E10 . 3 , 

10 _  218,*  ITERATIONS*//) 

11  FORMAT (1X,*A(SAMP.E) =* , 1 JE9. 3 , 3X, *K (S AMPLE) = * , £9.  I,5X,*N0=*,E10.3," 
218,*  ITERATIONS*//) 

_ 12__  FORMAT (2X, •INCOMPLETE, RE1JE5T  MORE  TIME*) 

i3  FORMAT  (  15X,*C0MPUTATIDN  I t ME=* , F6. 1 , •  SEC 3N05* / 2 IX , *PRE C ( A , K ) »* 
25  2.F8.5/I 

_ 14 _  FORMAT (19X, ’DILUTION2* , 1 3 El 0  •  3F) 

"16  FORMAT  (lH*,36X,*a£*0W  THIS  DIA,  N  IS"  UNCERTAIN.  IHCREASE  IN  P  IS  N 
2EGLIGI8LE.*) 

_ 17 _  F0RMAT(1H+,3oX,*AB0VE  THIS  JIA,  N  AND  P  APPRJADH  ZERO  RAPIOLT.*) 

'30  19  FORMAT  <//2X,*N0  SIGNIFICANT  ’ARTICLE  COUNT  H  THIS  SAMPLE,  */2X,*DL 

20W=*,I3,  3X ,*  OMIGH=* ,  I3/2X, *(SEE  TABLE  FOR  Z  2ALJSS  PRINTED  ON  NEXT 
3  PAGE)*) 

21  FORMAT  (*  N  ANO  P  *ELL"  RA* IOlT  ,  *", 1 3 , •  ITEMS  WiRE"S  ElITE  D  IN  TABlE*) 
1000  REA0(5,1)TITLE 

35  IF (EOF (5) ,NE.O.) STOP"NORMAl  EXECUTION” 

0  I L  5 1 . 

OMEGA=7.5 

I8LAUK=IZTAB*0 _ _ _ _ _ 

PREC= 1 . E-4 
N  =  5 

T  =  90,  _ _ _ _ 

RE AO ( 5 , OAT  A ) 

DO  102  1  =  1,  N 

READ(5,2> IV(I,J)f J»l,16)  _  _  _ 

S*0 

DO  100  J  =  .  15 

S  =  S*V(1,J)  __  _ _ _ _  _ _ _  _ _ 

IF(S.tQ.V(I,16))50  TO  102 
DO  101  J  =  1 ,16 

IV(J)=V(I,J)  _ _ _ _ _ _  _ _ 

HRITE(6,1)TITLE 
HRITE(6,3)  tlV(J), J*l,16) 

ST OP"CHECK  DATA” _ _ _ _ _  _ 

CONTINUE 
TIME=SECONO(S> 

CALL  OILSUB _ _  _  _ _ _  _ _ 

WRITE (6,1) TITLE 


40 


«*5 

100 

50 

101 

102 

55 
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PROGRAM  OIL 


TIM£  =  SFCONlHS>  -TIRE 
_  HRITE  <6,13)TIM£,3R£C 

IF(OM£GA.Gr. 0.150  TO  ZOOi 
IF  <  OMEGA.  EiJ.  0.  >HRm<6  ,H 
IF (OMEGA.LT.O.) WRITE (6,0) 

‘  STOP"£,E£  MESSAGE'* 

2001  IFUBLANK.tQ.OIGO  TO  3009 

_ _WRITE(b,10)  A2,F<2,FN0,IT£R___ _ 

IF<HOENO.E(I.U)Gd  TO  loOO 
WRITE (6,12) 

STOP~R£OUtST  MORE  T I  HE  “ 

3000  WRITE(6,14)0IL 

WRITE (6,11) A2,F<2, FRO, ITER 

IFINOENO. £0.0)30  TO  <*0  0  0  _ _ _ 

WRITE(6,12) 

STOP"REQU£ST  MORE  TIME" 

4000  IF ( IM, GT , 0 ) GO  TO  5000 

WRITE<6,19> IQLOdtORIGH 
IF (HI.5T ,0) RRITE(6,21I  RI 

_  _  _ CALL  ZTA3L-  _ 

GO  TO  1000 
5000  WRITE ( 6 , 4 ) 

_ 00  103  1*1, IM _ _ 

WRITE(6,5)iO(I)  ,D£LN(U  ,5EL»(I> 
IF(I.EQ.1)WRITE(6,16) 

103  __IF(I.tQ.IH.AN OiiM . LT.21MRITEI6, 1 7) 
WRITE ( 6, 6) 

IF(IZT4B.E0.1)CAL_  ZTABLE 
GO  TO  1000 


SUBROUTINE  0ILSJ8 
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subroutine  oilsub 
dimension  W  ( 15  1 

COMMON/lNPur/OI.,N,T,OME;A,V( 10, 16) ,IBLAN<,PREC 
C0MM0N/0UTPUT/UiLN(21)  , DEL3 ( 21 1 ,A2,FK2,FN0, If £<,11,13(211  , MI 
COMMON/SOLVfi/A3,F<8,FN03,Z, I JLOK , IOHIGH,<RI T 
COMMON/GRIOB/lTERArNOENT'  '  — 

U0MH0N/0CAL/0(ls,n 

COMMON/A UX/M,X (15)  ,  Y (16) , YO , DD, SUNY,FX (15> , YJ AL (l S) 

EXTERNAL  SUMDEO, SOLVER 

IM=MI=0 

DO  100  J -1 , 15 

00  100  1=1, N 

100  _ W(J)=W(J)+V(I,J^ _ _ 

V ( 16) =0 

00  101  1  =  1,15 

M=16-I 

Y(M)=Y(Moi)OH(H) 

101  CONTINUE 

_ M=  1  _ _ .  _ 

00  lu2  1=10,35,5 
M=M*1 

102  IF (OMEGA. EOi.FLOAT (I) )NQH«M 

if ( omega. eq.?.5>  n£)R=T  /  “ 

SUMY=M=I M=0 
DO  103  1=1,15 

If  (Y  ( 1 1 •  td, o .T&O  ro"  i'S'J 

M=M*1 

_  X(H)=0(M,N0M) _ 

Y(M)=AL0G(Y(N1) 

SUHY=SUMY+Y (N) 

103  CONTINUE  _ _ 

IF (M.Gt • 2) GO  TO  120 

0MEGA=-1, 

RETURN 

120  00=0(16, NON) 

ASTP=FKSTP=,03 

______  A  =  ,  25 _ 

FK=  1, 

PRECI=PREC 

_ _ _ CALL  GRI02 (SUHOEV, PRECI . 1ST » , FKSTP. A, FK.S) _ 

lF(S.NE.i,E99)50  T5  150 
OMEGA=0 

_ RETURN _ _ _ 

150  I TER=I TERA 
A2  =  A 
FK2=FK 

>N0  =  i.622*EXP(Y0i  /(FLOAMN)  *T> 

IF ( IBL  ANK. EQ, 0 1  GO  TO  121 

_ FNOB=FNO  _ _ _ _ _ 

AB  =  A 
FKB=FK 

_ I F  ( I  EH.  ANK,  EQil  7  RETURN _ 

121  IFINOENO.EQ.DRETJRN 

.  0ELNX0=-1.E99 

_ 00  200  1  =  1.95 _ _ _ _ _ 

IF ( I »GT, 11 OELNXO* OELNX 


35 


SUBROUTINE  OILSJB 


t  DELNX  =  FNO  *c.X®  (- A 2» FLOAT  ( I)  * •  FK2> 

2  -FNOU»£X®<-AB»FLOAT(  n ••FKS) 

>0 

200 

IF(0£LNX.LT.O.)GO  TO  200 

IF (DELNX. LT.UELNX0)G0  TO  201 

CONTINUE 

6$ 

201 

I0L0W=I-1 

lF{IOLON.Gi.94)  RETUIN 

ID(1)=IDLDH 

9000 

DL  OH: I DLQW 

or=oto^ 

OYsOY  +  U 

TO 

IF(JY.GT.95.)G0  To  999 

ZZs30LYER(0f|DUiY) 

IFIZZ.GT.Q.IGO  TO  9001 

9001 

OY=QY-l. 

GO  TO  999 

IF(ZZ.GI..675.0T.<?IT.EQ.1>GO  TO  90)0 

TS 

999 

luHIGH=QY 

IF  (IDLOH.EQ.IOHIGY.OR.  1 0-1 1 91.  LT.  5)  RETURN 

IFUOHlnH.GT.80)  IOH(GH*80 

80 

DHIGHUDHIGH 

DO  300  1=1,2 

IF(I.EQ.l)FLOATI  =  J*.QN 

300 

IF  (I.LO.2)FLOATI20HIG-( 

DELN(II*(FNO  »iXO(-A2»FtOATt  »*FK2>  ' 

2  -FN0B*EXP(-AB»Fl3ATI  **FKB)>*OIL 

85 

OSAVE=DElN(2> 

XHIB=AB»0HIGH*»F<3 

XMIS*A2*OHIGH»'FK2 

ARB=3. /FK9*1. 

ARS=3./FK2+1, 

PIb=ALOG(5.235983E“7) 

90 

C 

CALL  MU&AM(XHI9,A<U,RHI3iIE() 

CALL  NOGAM(XHIS,A<S,RHlS,IER> 

FOR  FUNCTION  NAMES  SEE  NEW  INSL  LIBRARY 

95 

FLOGHB=AL GAMA (ARB) 

FLOGMS=ALGANA»«RS) 

XLOH^AB'DLOH'yFO 

CALL  MOGAM<XlOW,ARB,RLOi(,IER> 

P1=0 

lF(kLOM«E3.-HI5)GO  TO  3900 

100 

Pl=P16+ALOG(FNoai  -3.*AL0J(A9)  /FK8tFl0GNB*AL05 ( ASS ( R.ON-RHI B> > 

IF(P1.LT,-8A. )Pl=-84 

IF(Pl.GT.Bb.)Pl  =  B» 

3000 

PB=E  XP (P 1) 

IF (PI.  LQ. 0. ) PB  =  0 
xlow=A2»olow**f<; 

105 

CALL  NOG  AN (XlOH|ARS,RLOM,IER> 

P2  =  0 

IFIRLOW.El.RHIBIGO  To  4000 

110 

P2=PI6+ALOG(FN0  ) -3. »AL0G<12> 7FK2»FL0GNS*AL0; (ASS ( R.ON-fcHIS) ) 

IF(P2.GT.8b. )P2=84 

IF(P2.LT.-84.)P2=-84 

9000 

PS=EXP(»2) 

IF(P2.ER.O.)PS*0 

0ELP(1)=0IL» <PS-?3) 

DEEP ( 2) >0 


SUBROUTINE  OILSJB 


t 

.  its  . 

*  ‘  y»»" 

oo  ini  i«i,z 

IFIOELNCD.LT. 0.»OELN<I)»0 

till  CONTINUE 

if  <  oel r.  i  r ;  ottit.  pc  rr*u 

IF  COELNC  I). El .0.. AN].  DElPC I > . EQ. 0 • >  1I*NI  *  1 

1Z0 

J=2 

00  108  1=5,95,5 
IF(I.LE.IULOW)GO  TO  108 

Tn'ircin  OHi  sh  rsiTTirs  5T 
j=j+i 


125 

DELNC J )  =  CFNO  »£XPC-A2»FL]ATC I) *»FKZ> 

* 

T  -FICUS*  Eli’T  -A  5*FlSITC  I  J**FKB) )  *DtL  ~ - - ~ 

IF  ( QELN  < J )  >.T,C,IOE.NCJI=0 

XL 0H  =  AS* FLOAT  CI)**FKB 

130 

CALL  HOG  AM (XLOH , A  98, 9L5iCT1?9> 

P1  =  0 

IF(SLOW.EQ.9HIB>SO.  TO  1000 

~~ Pi=PI6*Al06<FNQB) -3,*fCKTlS)  /FKB*FLtJCTB*ALOJ  USSTTLS > 

IFIPl.LT. -84. )P1=-8A 

IF (Pl.GT. 8b.) Pl=8 5 

135 

PB=EXP  (PI) 

1000  IF(P1.EQ.O.IPB=0 

XLOW=AZ*FLOAT<I> **fK2 _ 

CALL  MOG  AM  C  XLOH » A95  ,~RLOH»  I ET) 
P2=  U 


140 

IFCRLOM.EI.RNISIGO  TO  2000 

PZ  =  Pi6*AL0UCFNU  )  -3«*AL3G(A2) /FK2+FLOGMS+ALOI (A9S(9..]«-RMiS)  > 

IF(P2.GT.86,1P2=85 

IFCP2.LT. -04. )P2=-84 

J 

PS=EXP  CP2) 

• 

145 

2000  IF(PZ.Ed.0.)PS=0 

DELPIJ)=OIL*<PS-PB> 


IF(oeLP(JI.Lr.Q.)  DE.PCJ)  =0 

IF«OtLN(J)  .El.  0..  AND.  J-LPCJ)  .  EQ.  0  , )  M  I  =  M£  »1 

IFCOELNCJ) .E3.D..AN3.Di.P(J).EQ.O.)GO  TI  500 


* 

150 

108 

500 

ID  < J»  =  1 

CONTINUE 

IM=J 

IF  ( IM.GE.  3)  50  TO  70J 

IM=-IM 

155 

PET  URN 

700 

JM=J-1 

DO  501  J=2, JM 

IO(J)=IO(J+l) 

160 

501 

OELN(J)  =  OELN(J.l) 

DELPCJ)=DELPCJ»1) 

ID  C IM) =1 DHIGH 

OEL  P  C IM) =0 

OELN(IH)  *OSA VE 

00  600  1=1. IN 

165 

600 

OELPCI)=OELPCI)*5.2J59»»E-7»OELNCIM>*CF.OATCIOlI5H))  »*3 
RETURN 

ENO 

[ 


37 


SUBROUTINE  ZTABl  E 


1 


5 


10 


15 


SUBROUTINE  ZTA91E 

_ C0HM0N/S0LVB/A9,F<9iFM>9,4,  t  PLOW,  IPHIGHtKRIT 

WRITEiBf  32 > 

32  FORMAT ( 1  HI ,1 UX , ♦ T  ABL E  F09  Z*> 

MRITE<5, JO) 

30  FORMAT  «'/9X,*U»,TXi»Z»//) 

31  FORMAT  <110 ,F9.3) 

39  FORMAT (1H»  ,22X ,  »NS»NB>  10  ) » ) _ 

NO  FORMAT 11H+  »22X| *T ilS  IS  THE  50HPUTEQ  OHIGM*) 

00  $0  I* I DLOM 1 95 
OY=I 

S=50LVc9(DY ,0UXY) 

WRITE  (5>31)I|2 

_ _ IF(  IKrUI.EQ.l)  AUT-1S.39) _ 

so  iFa.ea.iOHisH)  orxte to, noi 

RETURN 

END 


FUNCTION  SOLVE R 


1 


5 


10 


FUNCTION  S01VER(0Y,8UMY1 

C0MM0N/0UTPUT/0£LN(2U  »DEL»(’l)  ,  A2,FKZ,FN0,IT  ER,I1»  13«21>  ,«I 
C0MM0N/S0LVB/A8f  F<B'» FNOBj  Z»  1 9LOH f  IOHIGH  jKRIT 
KRIT  =  0 

IF  (OY. LT. FLOAT (IOL ONI.  09. Of. JT. 95.) RETURN 
'FNS=FN0»  tXP<-A2»OY»»FK2T 
FNB=FN08*EXP<-Ad*3Y**FK3t 

IF<FNS+FNS.3T.10l*«>  <9IT»1 
l X  < FNS-FN3)  /SORT  ( FN5>FN 9 ) 

Z=FL0ATUFIX(Z*l.E3+,5))*t.E-3 

$OLVES=Z 


RETURN 

END 


FUNCTION  SUHDEV 


H 


r. 

& 


4 

— v - 

1 

FUNCTION  SUMOEtf ( A , FK1 

COMMON/AUX/N,X  (15)  ,7(161  ,  YO  ,  DO  ,SUNT,  FX(  161  ,Y:AL(161 

5 

SUM0EV*1.E99 

IF  (A.LE.  Q..OR.FK.LE.O.  1  RET  JIN 

SUMX=0 

00  100  1*1, M 

FX<Il=EXP(-A*X(i) **FK1  -£XP(-A»DD**FKl 

IF(FX(I).LE.D.1R£  TURN 

10 

FX(I1*AL0G(FX(II> 

100  SUHX=SUMX4FX ( II 

Y0=(SUMY-SUHX1 /F. DAT (HI 

SUHOEV*0 

00  101  1*1, M 

VCAL(I1*Y04FX(I) 

15 

101  SUHOEV*SUMOEYMY(Il-YCAL(II  1**2 

SUMOEV*FLOAT(IFIX (S JMO- V*1 . E9*. 51 1 *1. E-5 

RETURN 

BLOCK  OAT*  BLKOIT. 


I  .  1  BLOCK  OATA 

)  '  . _ C0MM0N/0CAL/0(16,71 

0ATA(0(l,Jl,  J*l,7i/2.07,2.15,2.34,2.53,2.74,2.  95,3.21/  “ 

i  UATA(U(2,Jl,Jsl,7)/5.22,5.4?,5.B6,6.32,b.83,7.37,  7,55/ 

_ _ 5  OATA  (0  (3,  J),  3*1,71  /10.67. 11. 0  7. 11. 94. 12.  86. 11. 87.  14.94,15.11/ 

*  DAlA(0(4,Jl, J*l, 71/15. 17, 16. 77, 14. 06,19.44,21. 93,  2  2,53,24.26/ ' 

DA TA(D (5, J1 ,J=1, 71 724. 36,25.26,27.16,29.20,31 .41, 33.76,36.32/ 

i  ..  __ _ _ OATA  (0  (6,  Jl,J*l, 7 1  732.  07,33,2  4,3  5.71,38.37,41.25,  44.3  2,47.52/ 

j  '  DATA (0 (7, j) ,J* 1,71  735.  26, 37. 57,40. 36, 43. 35, *5. 58, >0.0  4,53.  75/ 

10  OATA  (0(8,  J)  ,  3*1,71/3  9.  82, 41.26, 44.31, 47. 58, >1.11, 54.89, >8.  94/ 

_ _ _ OATA  (0(9, J),J*  1,71/43.01, 44.56, 47.94, 51.36, 55. 17, 59.23*63.59/ 

OAf*(0(id,Jl ,j*l, 73 745. 38, 47. 5 2, 51. 02, 54. 76, 38. 81, 53. 14,57.77/ 
OATA  (0  (11,  J),J»1, 7 1  748.  7  2,  53. 46, 54. 16, 58.  13, 5  2.  42, 65.  99, 7 1.90/ 

' - . - ~ - - - -  .  0ATA(0(12.J1  ,J»lL71j751_^37f  53.20,57.10,61.27, 55.78,70.  >9,75.75/ 

J  15  OATA ( D ( 1 3, Jl ,J*1, 71  754. 15, 5b, 08, 60. 17, 64, 56, 3 9. 30, 74.  37,79.79/ 

|  DATA (0 (14, J!,J»1, 71/5  5. 58, 5  6.  69, 62. 97,67. >6,7 2. >1,77.  80,83.46/ 

L_ _  DATA  (0  (15,  J)  ,J»1, 71  7  59.  43, 61. 56,66.04.70.84.76.02.81.55.47.48/ 

DATA ( D ( 16 , j>,J*l,7)/62,27,66.  47,69. lfe,74.  1 7, 79.53,8 5.37 ,51.577 
END 


SUBROUTINE  GRIDS 


• 

1 

SUBROUTINE  GRI02 C SUNCHI , if E»ll , XSTEP, YSTEP.X, Y,*l 
COMMON/GRIU8/ITiAA,NOEN3 

5 

NOENU  =  K=ITERA  =  tl 

PI =ATAN2  (0.,-l.) 

PI2=ATAN2(1.,0.) 

■ 

F  =  F0  =  FUNCMUX,Y> 

TIME  =  THETAO=l.t99 

1000  IFITIMREMU.) . LT. SECON) ( 1. > -TIKE*2. » NOEXO*t 

10 

T IHc=SECON3( l. ) 

1J=I1=J1=1 

ITtRA*  IT£RA*1 

IF  (XSTEP.LT*  ST  t!*llli  1*3 

IFIYSTEP.LT.STEPLI)  Jl*3 

IF  (N0LN0.EQ.1.0R.  <«E0.  20)I1*J1*3 

15 

IF (I1RJ1.LT. 6) GO  TO  100 

F=FUNCMI (X,Y) 

RETURN 

105  Xc=x 

YC  =  Y 

20 

FC=F 

DO  101  1=11,3 

FI.(C  =  FL0»TU-2>»X5TE» 

IF(n.LQ.3)FINC  =  0 

XG  =  XCFINC  / 

25 

DO  101  J- Jl,3 

FJNC=FLOAT(J-2)*Y5TEP 

IF(J1.EQ.3)FJMC=0 

IF (FINC. LQ. 0. • ANQ. F JNC. E). 0. ) GO  TO  101 

YG*YCfFJNC 

30 

F-G=FUNCMI  (XG.YG) 

IF (FG.GT.F1GO  TO  101 

F*FG  .  . 

X  =  XG 

Y  =  YG 

35 

THET  A=A1AN2(F  JNC , FINC) 

IJ=IJtl 

101  CONTINUE 

IF(F.LT.FC»GO  TO  102 

GO 

IF  (F.  EQ<  FO  <ANO. I J*3  * IAB5  (I1*J1).£Q»9)G0  TO  1)3 
F0=-1.E99 

K*0 

F*FC 

X*XC 

Y*YC 

G5 

XSTEP*. 5*XSTEP 

YSTEP=.5*YSTEP 

GO  TO  1000 

102  IF(THETA.EQ.THETAD)GO  T)  19J 

* 

* 

50 

K*fl 

_  THt T40=THtTA 

GO  TO  1000 

103  K=K»1 

IF (ABS  (THETA) .N^.PI2)XST-P*S.»XSTEP 

55 

IF  (THETA  .  NE.  PI »  0X0.  THE  TA.NE.O.)  TSTEP*2,  *f STE’ 

GO  TO  1000 

ENO 

40 


~«r*-  r-rn* 


jTofcfr  i  M0L  k  Z^.  j 


APPENDIX  B 

DATA  BASE  FOR  AUTHENTIC  OIL  SAMPLES 
Code  for  column  headings: 

T  =  operating  hours  elapsed  since  last  oil  change 

X  =  microliters  of  particulates  per  liter  of  oil, 
measured  by  particle  counting 

Y  =  microliters  of  particulates  per  liter  of  oil, 
inferred  from  spectrometric  analysis 


SAMPLE 

NUMBER  _ 

FAI LURES ( 31  ) 


ENGINE 

TYPE 


F-002 

F-100 

F-l  5 

219 

4.15 

5.81 

0.71 

F-003 

T-56-A7B 

C-130E 

113 

5.75 

2.35 

2.4 

F-004 

T-56-A7B 

C-130E 

493 

7.02 

5.56 

1.3 

F-005 

T-56-A7B 

C-130E 

755 

0.055 

11  .0 

<0.1 

F-006 

T-76-10A 

0V-10 

19 

14.4 

2.57 

5.6 

F-008 

T-56-A7B 

C-130E 

689 

90.7 

10.8 

8.4 

F-009 

J-69-25 

T-37B 

2 

25.8 

6.24 

4.1 

F-012 

T-56-A7B 

C-130B 

76 

14.5 

9.27 

1.6 

F-013 

TF-41-1 

A-70 

68 

6.85 

4.00 

1  .7 

F-  01 4 

T-56-A7B 

C-130E 

535 

1 .66 

8.67 

0.19 

F-015 

T-56-A7B 

C-130E 

unk 

158. 

30.9 

5.1 

F-017 

J-85-5 

T-38 

3 

26.6 

4.42 

6.0 

F-018 

T-56-A7B 

C-130E 

527 

125. 

35.9 

3.5 

F-020 

TF-33 

C-141A 

783 

10.6 

1.30 

8.2 

F-023 

T-56-A7B 

C-130E 

400 

307. 

31  .3 

9.8 

F-024 

T-56-A7B 

C-130E 

439 

0.387 

11.3 

<0.1 

F-025 

TF-33P7 

C-141 

814 

23.7 

3.78 

6.3 

F-026 

J-79  15 

RF-4C 

unk 

35.5 

29.4 

1.2 

F-028 

T-56-A7B 

C-130 

625 

17.5 

11.3 

1  .5 

F-030 

T-56-A7B 

C-130E 

790 

8.25 

5.76 

1  .4 

F-  031 

T-56-A7B 

C-130E 

52 

2.96 

6.75 

0.44 

F-033 

J-79-GE-1 5 

RF-4C 

51 

28.1 

6.57 

4.3 

F-034 

T-56-A7B 

C-130B 

284 

14.1 

6.81 

2.1 

F-035 

J-85-5 

T-38 

66 

6.83 

5.83 

1  .2 

F-036 

J-85-GE-5 

T-38A 

100 

29.5 

2.18 

14. 

F-037 

TF-39-GE-5 

C-5A 

617 

19.2 

11.6 

1  .6 

F-039 

T-56-A7B 

C-130E 

130 

26.1 

23.3 

1.1 

F-040 

T-56-A7B 

C-130E 

unk 

26.2 

10.1 

2.6 

F-041 

T-56-A7B 

C-130E 

1081 

74.6 

11.1 

6.7 

F-042 

T-56-A7B 

C-130 

1301 

41.5 

2.80 

15. 

F-046 

F-100-PW-100 

F-15A 

149 

17.7 

2.21 

8.0 

SAMPLE 

NUMBER 

ENGINE 

TYPE 

AIRC 

TYPE 

H I TS ( 24 ) 

H-002 

RT-6A 

C-130T 

H-009 

J-75-17 

F-106 

H-015 

J- 79-15 

F-4C 

H- 01  8 

J- 85-5 

T-38 

H-025 

J-57-59W 

KC-1  35 

H-044 

J-79-15 

F-4C 

H-048 

J-85-5 

T-38 

H-050 

T-56-A7B 

C-130 

H-052 

TF-39 

C-5A 

H- 06 1 

T-56-9 

AC-130A 

H-063 

J-85-GE-5 

T-38A 

H-065 

T-56-A7B 

C-130E 

H-070 

J-79-GE-1  5 

F-4C 

H-074 

J-85-GE-5 

T-38 

H-075 

J-85-GE-5 

T-38 

H-077 

J- 79-GE-l 5 

F-4C 

H-078 

TF-33-P-7 

C-141 

0-002 

J- 79-GE-l 5 

RF-4C 

P-002 

GE-T-58-5 

CH-3E 

P-025 

I0-360D 

0-2A 

P-050 

R-2800-1 03 

C-131D 

P-069 

T-400-CP-400 

UH-l-NA 

P-090 

J-85-GE-1 5 

T-38A 

P-091 

J-85-GE-15 

T-38  A 

ROUTINES (145) 

R-001 

TF-41 

A-7D 

R-002 

TF-41 

A-7D 

R-003 

T-56-A7B 

C-130E 

R-005 

T-56-A7B 

C-l  30 E 

R-007 

TF-33 

C-141 

R-011 

T-56-A7B 

C-130 

R-017 

T-56-7 

C-130 

R-024 

F-100DFCA 

F-15 

R-30 

J-79-15 

RF-4C 

R-033 

J-79-15 

RF-4C 

R- 041 

J-85-GE-5 

T-38A 

R-044 

J-57-59 

KC-1 35 

R-053 

J-79-15 

RF-4C 

R-056 

J-57-59W 

KC-135A 

R-058 

T-56-A7B 

C-130E 

R-077 

J-57-43W 

B-52G 

R-082 

J-57-59 

KC-135A 

R-084 

J-57-59 

KC-1 35 

R-086 

TF-41-A-1 

A- 70 

R-089 

J-79-15 

RF-4C 

R- 103 

J-57-43W 

B-52G 

T 

X 

Y 

X/Y 

18 

16.0 

15.4 

1  .0 

107 

25.2 

9.69 

2.6 

210 

4.57 

15.1 

0.30 

309 

21  .0 

4.61 

4.6 

134 

30.6 

9.27 

3.3 

84 

46.5 

14.0 

3.3 

unk 

110. 

163. 

0.67 

868 

10.6 

15.5 

0.69 

508 

34.8 

40.3 

0.86 

442 

64.9 

29.9 

2.2 

no 

4.19 

9.69 

0.43 

568 

15.8 

25.3 

0.62 

150 

1  .01 

10.0 

0.10 

41 

17.5 

4.52 

3.9 

138 

4.55 

6.72 

0.68 

550 

1.75 

9.52 

0.18 

576 

8.52 

10.5 

0.81 

221 

83.2 

17.0 

4.9 

51 

82.2 

49.8 

1.7 

44 

1270. 

77.7 

16. 

unk 

968. 

38.9 

25. 

99 

9.71 

9.25 

1  .1 

unk 

2.66 

9.95 

0.27 

unk 

10.9 

7.27 

1  .5 

56 

nil 

0.43 

<0.1 

24 

nil 

0.76 

<0.1 

272 

89.8 

20.5 

4.4 

310 

nil 

13.9 

<0.1 

1 

nil 

1.74 

<0.1 

156 

ni  1 

2.70 

<0.1 

136 

nil 

3.48 

<0.1 

15 

nil 

0.50 

<0.1 

450 

nil 

2.24 

<0.1 

69 

nil 

1  .79 

<0.1 

267 

nil 

3.49 

<0.1 

77 

nil 

2.36 

<0.1 

556 

nil 

2.15 

<0.1 

unk 

nil 

1.15 

<0.1 

59 

nil 

0.63 

<0.1 

179 

83.1 

2.67 

31  . 

0 

nil 

1 .45 

<0.1 

60 

nil 

2.75 

<0.1 

49 

nil 

6.91 

<0.1 

304 

nil 

2.04 

<0.1 

0 

nil 

4.01 

<0.1 

SAMPLE 

NUMBER 

ENGINE 

TYPE 

R- 104 

J-57-43W 

R-106 

TF-41 

R-122 

J-57-43wb 

R- 127 

J-57-43 

R- 129 

J-79-GE-15 

R-135 

J-57-59w 

R-148 

J-85-5 

R-149 

J-85-5 

R- 153 

J-57-43wb 

R-164 

J-57-59w 

R-155 

J- 79-1 5 

R- 157 

J-57-43W 

R-158 

J-57-59 

R- 161 

J- 79-15 

R-164 

J-57-59 

R-167 

J-79-15 

R-177 

J- 79-1 5 

R-178 

T - 56- 7 

R-183 

J-57-43wb 

R-186 

J-57-59W 

R-190 

J-69-T-25 

R- 199 

J-57-59W 

R-200 

J-57-59w 

R-201 

J-57-59 

R-203 

J-57-59 

R-204 

J-57-59 

R-206 

J-79-15 

R-215 

TF-41 -1 

R-219 

J-57-43wb 

R-220 

J-57-59w 

R-226 

J-57-59w 

R-227 

T-56-7 

R-229 

J-57-43wb 

R-238 

J-57-59w 

R239 

J-57-43W 

R-240 

J-79-15 

R-241 

J-57-59w 

R-245 

J-57-59W 

R-246 

J-57-43wb 

R-247 

J-57-43wb 

R-273 

J-79-15 

R-274 

T-56-7 

R-283 

J-57-43wb 

R-298 

J-85-5 

R-299 

J-69-25 

R-300 

J-85-5 

R-311 

J-57-59W 

R-315 

J-57-43wb 

R-316 

0-57-59W 

AIRC 

TYPE 

T 

X 

KC-135 

1132 

ni  1 

A- 70 

71 

nil 

B- 52 

42 

nil 

B-52G 

220 

nil 

RF-4C 

127 

nil 

KC-135 

6 

nil 

T-38A 

22 

nil 

T-38A 

144 

nil 

B-52G 

1652 

nil 

KC-135 

703 

nil 

RF-4C 

16 

nil 

B-52G 

408 

nil 

KC-135 

900 

nil 

RF-4C 

571 

nil 

KC-135 

141 

nil 

RF-4C 

16 

nil 

RF-4C 

661 

nil 

C-130 

27 

nil 

8-52G 

812 

nil 

KC-135 

300 

ni  1 

T-37B 

21 

nil 

KC-135 

593 

nil 

KC-135A 

107 

nil 

KC- 1 35  A 

505 

nil 

KC-135A 

328 

nil 

KC-135A 

291 

nil 

RF-4C 

412 

nil 

A-7D 

92 

nil 

B-52G 

300 

nil 

KC-135 

10 

nil 

KC-135A 

462 

1.77 

C-130 

380 

nil 

B-52G 

1169 

nil 

KC- 1 35 A 

125 

nil 

B-52G 

277 

nil 

RF-4C 

6.3 

ni  1 

KC-135 

9 

nil 

KC-135 

526 

ni  1 

B-52G 

1434 

nil 

B-52G 

893 

nil 

RF-4C 

163 

nil 

C-130 

61 

nil 

B-52G 

300 

nil 

T-38A 

250 

8.31 

T-37 

21 

ni  1 

T-  38 

1 

nil 

KC-135 

869 

36.0 

B-52G 

304 

nil 

KC-135 

102 

nil 

_Y _ 

2.00 

3.66 

4.44 
2.03 
2.24 
3.33 
2.47 
2.88 
2.30 
2.64 
0.43 
1.39 
1  .90 
1  .92 
1  .27 
0.91 
2.00 
2.35 
3.11 
1  .22 
0.87 
2.01 
2.18 
4.26 
1.28 
2.24 
2.04 
2.96 
0.45 
1 .33 
12.1 
2.65 

1 .45 
1 .82 
1  .34 
2.31 
0.71 
3.68 
5.91 
3.73 

1.25 

4.45 
0.78 
8.93 
1.30 

1.26 
2.59 
1  .96 
3.03 


<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

0.15 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

<0.1 

0.93 

<0.1 

<0.1 

14. 

<0.1 

<0.1 
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SAMPLE 

NUMBER 

ENGINE 

TYPE 

AIRC 

TYPE 

T 

X 

Y 

X/Y 

R-317 

0-57-43wb 

B-53G 

300 

nil 

1 .03 

<0.1 

R-  31 8 

J-85-5 

T-38A 

102 

7.70 

7.56 

1.0 

R-319 

J-79-15 

RF-4C 

318 

nil 

2.75 

<0.1 

R-324 

T- 56- 7 

C-130 

1166 

5.27 

11.4 

0.46 

R-325 

T-56-A7B 

C-130E 

11 

3.33 

4.07 

0.82 

R-340 

J-57-P-43wb 

B-53G 

985 

nil 

1.72 

<0.1 

R-342 

0-79-15 

RF-4C 

12 

nil 

1 .14 

<0.1 

R-345 

J-57-P-43w 

B-52 

222 

nil 

1.65 

<0.1 

R-  353 

J-85-GE-21 

F-5E 

3 

nil 

0.66 

<0.1 

R-366 

0-85-GE-5 

T-38A 

91 

nil 

2.45 

<0.1 

R-372 

J-79-GE-15 

RF-4C 

107 

nil 

1.45 

<0.1 

R-386 

J-69-T-25 

T-37B 

80 

nil 

1 .85 

<0.1 

R-387 

0-85-GE-5 

T-38A 

42 

nil 

1.99 

<0.1 

R-390 

J-85-GE-5 

T-38A 

109 

nil 

2.09 

<0.1 

R-  391 

J-69-T-5 

T-37B 

168 

nil 

1 .56 

<0.1 

R-400 

J-57-P-43wb 

B-52G 

972 

nil 

2.89 

<0.1 

R-404 

J-57-P-43w 

B-52G 

149 

nil 

3.96 

<0.1 

R-423 

J-57-P-59w 

KC-135 

615 

nil 

3.98 

<0.1 

R-424 

T-56-A7B 

C-130E 

2034 

49.6 

16.0 

3.1 

R-429 

J-85-GE-5 

T-38A 

150 

5.36 

5.11 

1.0 

R-433 

J-79-GE-15 

RF-4C 

572 

0.767 

14.3 

<0.1 

R-435 

J-57-P43W 

B-52G 

0 

nil 

1 .67 

<0.1 

R-436 

J-57-P-43W 

B-52G 

0 

nil 

2.12 

<0.1 

R-437 

J-52-P-59W 

KC-135A 

817 

nil 

2.60 

<0.1 

R-438 

0-57-P-43W 

B-52G 

1179 

nil 

1.78 

<0.1 

R-446 

J-69-T-25 

T-37B 

1 

nil 

4.77 

<0.1 

R-452 

J-85-GE-15 

T-38A 

unk 

348. 

3.14 

no. 

R-454 

J-57-43wb 

B-52G 

90 

1 .20 

3.62 

0.33 

R-455 

J-57-43wb 

B-52G 

2731 

2.39 

6.11 

0.39 

R-456 

J-57-43wb 

B52G 

32 

11.5 

2.64 

4.4 

R-457 

J-57-59w 

KC-135 

unk 

3.73 

7.83 

0.48 

R-460 

J-57-P-59W 

KC-135 

3405 

12.0 

3.03 

4.0 

R-461 

J-57-P-59W 

B-52G 

1613 

0.911 

3.02 

0.30 

R-462 

J-57-P-43 

B-52G 

1039 

nil 

1.6. 

<0.1 

R-465 

J-57-43wb 

B-52G 

3099 

nil 

2.53 

<0.1 

R-466 

J-57-43wb 

B-42G 

943 

nil 

2.74 

<0.1 

R-468 

J-57-P-59 

KC-135A 

687 

2.61 

3.26 

0.80 

R-473 

J-57-P-59W 

KC-135 

82 

3.21 

2.81 

1.1 

R-489 

T-56-A7B 

C-130E 

39 

nil 

3.88 

<0.1 

R-497 

J-57-P-59-W 

KC-135 

776 

5.97 

4.39 

1.4 

R-500 

J-85-GE-5 

T-38A 

192 

20.2 

10.8 

1  .9 

R-501 

J- 79- GE- 15 

RF-4C 

183 

2.67 

4.12 

0.65 

R-502 

J-79-GE-15 

RF-4C 

257 

1.21 

4.28 

0.28 

R-503 

J-85-GE-5 

T-38 

6 

25.4 

10.6 

2.4 

R-504 

J-85-GE-5 

T-38A 

136 

0.260 

8.14 

<0.1 

R-505 

J-85-GE-5 

T-38A 

unk 

7.04 

2.85 

2.5 

R-506 

J-57-43wb 

B-52G 

99 

0.922 

1 .76 

0.52 

R- 52  3 

J-57-P-59w 

KC-135 

unk 

nil 

1.53 

<0.1 

R-524 

J-57-P-59w 

KC-135 

313 

nil 

1 .90 

<0.1 

R-525 

J-85-GE-5 

T-38A 

210 

29.0 

6.92 

4.2 

R-527 

J-85-GE-5 

T-38A 

274 

1.03 

7.84 

0.13 

R-531 

0-79-GE-15 

F-4C 

437 

nil 

2.69 

<0.1 
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SAMPLE 

NUMBER 

ENGINE 

TYPE 

AIRC 

TYPE 

T 

X 

Y 

X/Y 

R-535 

J-57-43wb 

B-52G 

268 

nil 

1.87 

<  0.1 

R-541 

J-57-P-59 

KC-135A 

396 

nil 

3.26 

<  0.1 

R-542 

J-69-T-25 

T-37B 

unk 

nil 

1 .84 

<  0.1 

R-550 

J-79-GE-1 5 

RF-4C 

281 

1.48 

4.00 

0.37 

R- 551 

J-85-GE-5 

T-38 

200 

0.628 

4.85 

0.13 

R-552 

J-57-59w 

KC-135 

878 

nil 

2.36 

<0.1 

R- 553 

J-57-P-59w 

KC-1  35 

1485 

1.42 

2.53 

0.56 

R-554 

J-85-GE-21 

F-5E 

198 

0.523 

4.55 

0.12 

R-555 

J-69-T-25 

T-37B 

unk 

0.894 

6.82 

0.13 

R-556 

J-57-P-59W 

KC-135 

1252 

4.74 

5.16 

0.92 

R-  561 

J-79-GE-15 

RF-4C 

152 

0.782 

8.76 

<0.1 

R-562 

J-69-T-25 

T-37B 

unk 

0.652 

4.14 

0.16 

R-563 

T-56-A7B 

C-130E 

696 

4.48 

13.3 

0.34 

R-564 

J-69-T-25 

T-37 

unk 

4.20 

5.85 

0.72 

R-565 

T-56-A7B 

C-130E 

unk 

14.4 

2.67 

5.4 

R-566 

T-69-T-25 

T-37 

unk 

3.64 

3.39 

1.1 

R-570 

J-69-T-25 

T-37B 

135 

3.88 

3.79 

1 .0 

R-572 

J-85-GE-5 

T-38A 

141 

nil 

1.77 

<0.1 

F-032 

T-56-A7B 

C-130E 

27 

2.09 

4.25 

0.49 

F-038 

T-56-A7B 

C-130E 

765 

2.77 

6.28 

0.44 

H-031 

J- 85-2 1 

F-5 

150 

5.80 

6.74 

0.86 

H-062 

J-79-GE-15 

F-4C 

179 

0.109 

11  .1 

<0.1 

H-064 

J-85-GE-5 

T-38A 

2 

1  .47 

3.16 

0.46 
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APPENDIX  C 


THIN  FILM  PARTICLE  COUNTING 

Here  we  analyze  the  feasibility  of  performing  particle  counts  on  a 
thin  layer  of  undiluted  oil.  Based  on  experience  in  counting  diluted  oil 
samples,  we  show  how  the  allowable  film  thickness  is  related  to  particle 
concentration,  opacity,  and  counting  time. 

Particle  concentration 

According  to  our  present  measurements,  the  most  heavily-charged  oils 

may  contain  as  many  as  100,000  particles  per  cm3  (larger  than  one  micron 

in  diameter).  In  the  conventional  mode  of  operation,  such  a  specimen  must 

be  diluted  to  about  70  times  its  original  volume  in  order  to  reduce  the 

number  density  of  particles  below  1500  cm'^;  we  recall  that  this  is  the 

tolerable  limit  for  a  five-centimeter  path  length  if  multiple  scattering 

is  to  be  avoided.  The  characteristic  parameter  governing  this  behavior 

is  the  number  of  particles  per  unit  of  cross-sectional  area,  and  the 

-2 

critical  value  is  seen  to  be  5  x  1500  =  7500  cm  .  It  follows  directly 
that  the  maximum  permissible  thickness  for  an  undiluted  heavily-charged 
oil  is  L  =  7500/100,000  =  0.075  cm  =  0.75  mm.  Since  few  oils  approach  a 
concentration  of  100,000  cm*3,  this  thickness  is  tolerable  for  any  real 
specimen;  for  most  samples,  it  is  well  within  the  limit. 

Opacity 

At  a  dilution  of  70,  the  opacity  of  even  the  darkest  oils  never  exceeds 
35%  in  a  five-centimeter  cell;  at  least,  no  case  of  a  higher  opacity  has  yet 
been  discovered,  and  the  typical  oil  runs  much  lower.  We  wish  to  estimate 
the  opacity  of  such  a  dark  oil  in  the  undiluted  state.  Applying  the 
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Beer-Lambert  law,  we  write 


I/I0  =  1-0.35  =  0.65  =  0.925  exp  (-kLc) 

for  the  diluted  sample.  The  factor  0.925  in  this  expression  is  the  transmission 
of  a  glass  or  quartz  cell  containing  pure  paraffin  oil,  and  represents  the 
loss  due  to  reflection  at  the  cell  walls;  the  paraffin  itself  absorbs 
practically  no  radiation  at  the  laser  wavelength.  The  constant  k  is  the 
absorption  coefficient  of  the  undiluted  engine  oil.  If  we  supply  the  values 
L  =  5  cm  and  c=  1/70,  we  find  k  =  4.94  cm“^.  Using  this  constant,  we 
calculate  the  opacity  of  a  layer  of  undiluted  oil  0.075  cm  thick  to  be  almost 
exactly  35%.  Thus  we  see  that  even  an  exceptionally  dark  specimen  will  not 
exceed  the  desirable  opacity  limit. 

Counting  time 

Finally,  we  compare  the  counting  rates  for  undi  jted  and  diluted  samples. 
The  thickness  of  the  undiluted  sample  is  assumed  to  0.075  cm.  Recall  that 
the  sensitive  zone  of  the  counter  is  about  2  cm  long,  and  this  zone  is  fully 
utilized  when  a  diluted  sample  is  counted.  If  the  dilution  is  d,  then  the 
ratio  of  counting  rates  (undiluted/diluted)  is  simply  0.075d/2  =  0.037d. 

Hence  for  the  typical  oil,  run  at  the  normal  dilution  of  about  70,  the  count 
rate  in  the  undiluted  mode  will  be  about  2.6  times  as  fast  as  when  diluted; 
in  other  words,  a  count  time  of  180  seconds  produces  as  many  scattering 
events  as  occur  in  450  seconds  in  the  diluted  oil. 

In  the  case  of  an  oil  sparsely  populated  wi  '  particles,  which  in  the 
ordinary  way  might  sometimes  be  run  at  a  dilution  of  35  or  so  (opacity 
permitting),  the  thin-film  counting  time  would  have  to  be  doubled  to  realize 
equivalence  in  the  total  number  of  counts.  An  alternative  choice  would  be 
to  use  a  thicker  layer  (0.15)cm  and  retain  the  3-minute  counting  time.  However, 


it  is  questionable  whether  either  of  these  refinements  would  in  fact  be 
warranted;  the  principal  motive  in  reducing  dilution  during  a  conventional 
run  is  to  improve  the  blank  correction,  particularly  at  large  particle 
diameters,  and  this  problem  does  not  arise  in  thin-film  measurements. 
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